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Abstract 

In this paper, we consider the peak-covariance stability of Kalman filtering subject to 
packet losses. The length of consecutive packet losses is governed by a time-homogeneous 
finite-state Markov chain. We establish a sufficient condition for peak-covariance stability 
and show that this stability check can be recast as a linear matrix inequality (LMI) feasibility 
problem. Comparing with the literature, the stability condition given in this paper is invari¬ 
ant with respect to similarity state transformations; moreover, our condition is proved to be 
less conservative than the existing results. Numerical examples are provided to demonstrate 
the effectiveness of our result. 
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1 Introduction 

Networked control systems are closed-loop systems, wherein sensors, controllers and actuators 
are interconnected through a communication network. In the last decade, advances of modern 
control, micro-electronics, wireless communication and networking technologies have given birth 
to a considerable number of networked control applications. 

In networked control systems, state estimation such as using a Kalman filter is necessary 
whenever precise measurement of the system state cannot be obtained. When a Kalman filter 
is running subject to intermittent observations, the stability of the estimation error is affected 
by not only the system dynamics but also by the statistics of the packet loss process. The 
stability of Kalman filtering with packet drops has been intensively studied in the literature. 
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In [IHS], an independently and identically distributed (i.i.d.) Bernoulli packet loss has been 
considered. Some other research works assume the packet drops are due to the Gilbert-Elliott 
channel [6l[7], which are governed by a time-homogeneous Markov chain. Huang and Dey [8] 
introduced the notion of peak covariance, which describes the upper envelope of the sequence of 
error covariance matrices for the case of an unstable scalar system. They focused on its stability 
with Markovian packet losses and gave a sufficient stability condition. The stability condition 
was further improved in [UlIO]. In m, the authors proved that the peak-covariance stability 
implies mean-square stability for general random packet drop processes, if the system matrix 
has no defective eigenvalues on the unit circle. In addition to the peak-covariance stability, the 
mean-square stability was considered for some classes of linear systems in min], and weak 
convergence of the estimation error covariance was studied in M- 

In the aforementioned packet loss models, the length of consecutive packet losses can be 
inhnitely large. In contrast, some works also consider bounded packet loss model, whereby the 
length of consecutive packet losses is restricted to be less than a finite integer. A real example of 
bounded packet losses is the WirelessHART (Wireless Highway Addressable Remote Transducer) 
protocol, which is the state-of-the-art wireless communication solution for process automation 
applications. In WirelessHART, there are two types of time slots: one is the dedicated time slot 
allocated to a specific held device for time-division multiple-access (TDMA) based transmission 
and the other is the shared time slot allowing contention-based transmission. A contiguous group 
of time slots during a constant period of time forms a superframe, within which every node is 
guaranteed at least one time slot for data communication. Various networked control problems 
with bounded packet loss model have been studied, e.g., |151II6j : while the stability of Kalman 
hltering with the bounded packet loss model was rarely discussed. In the authors gave 
a hrst attempt to the stability issue related to the Kalman hltering with bounded Markovian 
losses. They provided a sufficient condition for peak-covariance stability, the stability notion 
studied in [EHio]. Their result has established a connection between peak-covariance stability and 
the dynamics of the underlying system and the probability transition matrix of the underlying 
packet-loss process. In this paper, we consider the same problem as in m and improve the 
condition thereof. The main contributions of this work are summarized as follows: 

1. We present a sufficient condition for peak-covariance stability of the Kalman hltering 
subjected to bounded Markovian packet losses. Different from that of this stability 
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check can be recast as a linear matrix inequality (LMI) feasibility problem. 

2. We compare the proposed condition with that of m- We show both theoretically and 
numerically that the proposed stability condition is invariant with respect to similarity 
state transformations, while the one given in m may generate opposite conclusions under 
different similarity transformations. Moreover, the analysis also suggests that our condition 
is less conservative than the former one. 

The remainder of the paper is organized as follows. Section [2] presents the mathematical models 
of the system and packet losses, and introduces the preliminaries of Kalman filtering. Section [3] 
provides the main results. Comparison with m and numerical examples are presented in 
Section 01 Some concluding remarks are drawn in the end. 

Notations: N is the set of positive integers and C is the set of of complex numbers. S” 
is the set of n by n positive semi-definite matrices over the field C. For a matrix X G 
cr{X) denotes the spectrum of X, i.e., (7{X) = {A : det(A/ — X) = 0}, and p{X) denotes the 
spectrum radius of X. ||X|| means the L 2 -norm on C” or the matrix norm induced by L 2 -norm. 
The symbol ® represents the Kronecker product operator of two matrices. For any matrices 
A, B, C with compatible dimensions, we have vec{ABC) = (C (8> A)vec{B), where vec(-) is 
the vectorization of a matrix. Moreover, the indicator function of a subset M C n is a function 
: Q ^ {0,1} where l^(w) = 1 if w € M, otherwise = 0. 


2 Problem Setup 

2.1 System Model 

Consider the following discrete-time LTI system: 

Xk+i = Axk + Wk, ( 1 ) 

Vk = Cxk + Ufc, (2) 

where A € and C G Xk G M” is the process state vector, yk G is the 

observation vector, Wk G M"" and Vk G are zero-mean Gaussian random vectors with 
K[wkWj'] = SkjQ {Q > 0), K[vkVj'] = 6kjR {R > 0), K[wkv/] = 0 Vj,k. Note that 6kj is 
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2.2 Bounded Markovian Packet-loss Process 


the Kronecker delta function with 6kj = 1 ii k = j and 5kj = 0 otherwise. The initial state xq 
is a zero-mean Gaussian random vector that is uncorrelated with Wk and Vk- Its covariance is 
So > 0. It can be seen that, by applying a similarity transformation, the unstable and stable 
modes of the LTI system can be decoupled. An open-loop prediction of the stable mode always 
has a bounded estimation error covariance, therefore, this mode does not play any key role in 
the problem considered below. Without loss of generality, all eigenvalues of A are assumed to 
have magnitudes not less than 1. We also assume that (A, C) is observable and (A, is 

controllable. We introduce the definition of the observability index of (A, G), which is taken 
from [T8] . 

Definition 1 The observability index Iq is defined as the smallest integer such that 

[C, A'C ,..., {A^°~^)'C']' has rank n. If Iq = 1, the system (A, C) is called one-step observable. 

2.2 Bounded Markovian Packet-loss Process 

In this paper, we consider the estimation scheme, where the raw measurements of the sensor 
{yk}k£N are transmitted to the estimator over an erasure communication channel: packets may 
be randomly dropped or successively received by the estimator. Denote by a random variable 
7 fc G {0,1} whether or not is received at time k. If jk = Ij h indicates that arrives 
error-free at the estimator; otherwise 'jk = 0. Whether jk equals 0 or 1 is assumed to have 
been known by the estimator before time A; -|- 1. In order to introduce the packet loss model, we 
further define a sequence of stopping times based on {7fc}fcgN) which presents the time instants 
at which packets are received by the estimator: 

ti = min{A: : k G N,'yk = 1}; 
t 2 — min{/c : A: > ti, 7 fc = 1}, 

; ( 3 ) 

tj = min{A: : k > kj^i^jk = 1}) (4) 

where we assume to = 0 by convention. Then at the jth. time instant the estimator successfully 
receives a measurement from the sensor. The packet-loss process, Tj, is defined as 

D ~ ^3 ~ 
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2.3 Kalman Filtering with Packet Losses 


As for the model of packet losses, we assume that the packet-loss process is modeled by 

a time-homogeneous ergodic Markov chain, where <S = {0,..., s} is the finite-state space of the 
Markov chain with s being the maximum length of consecutive lost packets allowed. Here this 
Markov chain is characterized by a known transition probability matrix 11 = in which 

TTij = P(Tfc+i = j\Tk+l =i)>0. (5) 

Denote by vr = [tto, ..., tt^] the stationary distribution of {rjjjgH- Without loss of generality, we 
assume that the initial distribution is P(ti = j) = ttj and other cases can be treated in the same 
manner. 

2.3 Kalman Filtering with Packet Losses 

Sinopoli et al. [T] shows that, when performed with intermittent observations, the optimal linear 
estimator is a modified Kalman filter. The modified Kalman filter is slightly different from the 
standard one in that only time update is performed when the data packet is lost at that time. 

Define the minimum mean-squared error estimate and the one-step prediction at the estimator 

respectively as Xk\k = ^[xk\liyi, ■ ■ ■ .IkVk] and Xk+i\k = ^[xk+ibwi, ■ ■ ■ .IkVk]- Let Pk\k and 
Pk+i\k Le the corresponding estimation and prediction error covariance matrices, i.e., 

Pk\k - ^[{xk-Xk\k){-)'\liVu-■ ■ ,lkVk] 

Pk+i\k = E[(xfc+i-Xfc+i|fc)(-)'|7i2/i,---,7fc2/fc]- 

These parameters can be computed recursively by a modified Kalman hlter (see [1] for more 
details). In particular, 

Pk+l\k = ^Pk\k-l^ + Q (6) 

-lkAPk\k-iC'{CPk\k-iC' + R)-^CPk\k-iA'. 

To simplify notations, we denote Pk — Pk\k-i shorthand and define the functions h, g, 
and g^-. §” —)• §” as follows: 

h{X) 4 AXA' + Q, ( 7 ) 

g{X) ^ AXA' + Q - AXC'{CXC' + R)-^CXA', ( 8 ) 

h^iX) = h o h o ■ ■ ■ o hiX) and g^iX) = g o g o ■ ■ ■ o giX), where o denotes the function com- 

'- 7 -^ '---^ 

k times k times 

position. 
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2.4 Problems of Interest 


2.4 Problems of Interest 

To study the stability of Kalman filtering with packet losses, one way is to study the asymptotic 
behavior of the expected prediction error covariance sequence. In the following we introduce 
the concept of peak-covariance stability, which is first studied in [8]. To this end, we need the 
following auxiliary definitions, which are also introduced in [5], 

= min{A; :/c G N, 7 fc = 0}, 

/3i = mm{k : k > ai,jk = I}, (9) 

aj = minj/c : k > /3j_i,7fc = 0}, 

/3j = min{A: : k > aj,jk = l}j 

where /3o = 0 by convention. It is straightforward to verify that 
sequence of stopping times (cf., [H]). 

Definition 2 The Kalman filtering system with packet losses is said to be peak-covariance stable 
z/supjgp,E||P^J| < oo. 

Note that E||Py 3 ^.+i||, the mean of one-step prediction error covariance at stopping time fij, reflects 
the stability of Kalman hltering at packet reception times. In the literature, stability of Kalman 
filtering with binary Markovian packet losses (driven by a two-state Gilbert-Elliott packet loss 
model) [8l[T0l[T2] and with i.i.d. packet losses [Tl[3] has been intensively studied. The main 
problem of this work is to study stability of Kalman filtering with bounded Markovian packet-loss 
process. As the packet loss is modelled differently, the stability also behaves differently. Due to 
the nonlinearity of the Kalman filter, it seems challenging to find necessary and sufficient stability 
conditions for a general LTI system. In Section [3j we manage to give a sufficient peak-covariance 
stability condition for general LTI systems with bounded Markovian packet-loss process. Our 
result is mainly built on the prior work HZ]. Compared with the result thereof, ours prevails 
from at least two aspects. We will discuss in details later in Section HI 


( 10 ) 


{ajIjgN and are two 
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3 Main Result 


In the following theorem, we will present a sufficient condition for peak-covariance stability of 
Kalman filtering with bounded Markovian packet-loss process. 

Theorem 1 Consider the system described in © and ([2]), and the bounded Markovian packet- 
loss process described by a probability transition matrix 11 m ([5]). If there exists K = ..., 

where ’s are matrices with compatible dimensions, such that p{Hk) < 1, where 

Hk = diag(^T T)") [P' ® H + Q' O K] , (11) 

and 

H ^ {A + KWC) ^{A + 

lo-l _ 

K ^ ^(7roo)^"^(A« + Kd)C(^)) ® (A^ + 

1=2 

P = Q — [^* 0 ^ 0 i]i,jG<s/{ 0 }/ then the state estimator is peak-covariance stable, 

i.e., supjgNE||P; 3 j| < cx). 


Before proceeding to the proof, we first present a few supporting definitions and lemmas. 
Consider k compositions of g together. We introduce the following lemma. 


Lemma 1 Consider the operator 

= {A^+ K^^C^^)X{-)*+ [A^^ 
\/i € N, 


QW QW(dW)' 




where C^^ = [C, A'C, • • • , {A'f-^C']', A® ^ [A^-\ ■■■ ,A, I], = Q for i = I otherwise 

0 0 0 


J)(i) A 


C 0 


0 


, = diag(Q,-- - ,Q), = diag(i?,-- - ,R), and 


CA^-^ CA^-^ ••• 0 

are of compatible dimensions. For any X > 0 and the following statement always holds 

g\X) = min<^,(l^«,X) < 
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Proof. The result is readily established when setting B = I in Lemmas 2 and 3 in m- For 
i = 1, The result is well known as Lemma 1 in [T]. □ 

The following lemma is about the nonlinearity of g operator: for k > lo+l, g^{X) is uniformly 
bounded no matter what the postive semidefine matrix X is. 

Lemma 2 (Lemma 5 in |8]) Assume that {A^C) is observable and is controllable. 

Define 

§0 = {T : 0 < P < APqA' + Q, for some Pq > 0}, 

Then there exists a constant L > 0 such that 


(i). for any X € Sg, g^{X) < LI for all k > Iq/ 
(a), for any X € g^^^{X) < LI for all k > Iq. 


According to the definitions of aj and flj, we can fnrther define the sojourn times at the 
state 1 and 0 respectively as follows 

a* A aj - /3j_i G N, 

P* = I3j - Oj G {1,... ,s}. 

The distribution of sojourn times aj and /3t is given in the following lemma. 


Lemma 3 (Lemma 4 in jl7j l Denote the joint distribution of a.*^ and /3t by 


tt{1) a p(ai = = bi,...,a^ = au(3f = k), 


for any Oj G N and fij G {1,... , s}. Then it holds that 


7r(l) 


7t{1 + 1 ) 


TT^l , if (^1 — 1; 

7ro(vroo)“^“^7rofei, ifai>2, 

^bibi^l'^{l)i if 0.1+1 — 1) 

i/a;+i>2, 



Proof of Theoreml^ We compute ]E[P^^] as follows: 


= E E 

ai=l 6 i=l 
s 

= o gi^o 

fei=i 


< 


lo s 

+ E E ’^o(7roo)“^“^7ro6i/l'’l 0 5r“l(So) 

0 , 1=2 6 i=l 
OO S 

+ E E ^o(^oo)“^"^7rofei/i^i o5“i(So) 

ai=lo+l bi=l 

s 

7rb,A’’^{A + K^^^C)J:o{A + K^^'^C)*{A^^y 

fei=i 


+ E' 4 ‘‘ 


fei=i 

s 


lo 

E vro(7roo)“i-Vofe,(^“i +i^(“i)C'(“i))So(-)* (A^^Y 
. 01=2 

Y ^b,A^^ K^^YJi[A K^Y) (A^^Y 

bi=i 

S lo 


bi=l oi =2 

s / \o \ bi-1 

+ E ( E ’^o(7roo)“^“^7robi + irbA Y A'QiA'Y 

bi=l \oi =2 / i =0 

OO s 

+ E E 7ro(vroo)“^ ‘^irobih’’^ {\-I) 

ai=lo+l fei=l 

= Ai + A2 + A3 + A4 + A5 + Ag, 


( 12 ) 


where J,- = 


and the inequality is from Lemmas [T] and [2l 


gW gW(nW)' 

DW(gW) D«(g«)(DW)' + i?« 

One can verify that A 3 , A 4 , A 5 and Ag are all bounded matrices. Then U = A 3 + A 4 + A 5 + Ag 
is also bounded. To facilitate the analysis in the following, we will impose ()12h to take equality. 
Without loss of generality, the conclusions in this paper still hold for other cases as (1121) renders 
us an upper bound of E[P^J. Next we vectorize both sides of (fT^ . One has 

s Iq 

E[vec(P^J] = A)^i [ ^ 7 ro( 7 roo)“i-Vobi(A“i + A:(“i)C(“i)) ® A) + 


bi=l 0 , 1=2 

7 rf,^(A + KA)C) (g) (A + kA')C) vec(Eo) + vec([/) 


(13) 


9 










Similarly, for any j > 1, E[P/ 3 ^_^J can be calculated as 

oo s oo s 

= EE'" E E r-ft+ii' + iwi + i) 

ai=lbi=l ajj^i=l bjj^i=l 
oo s oo s 

= + 1 ) 
ai=lbi=l aj+i=lo 6j+i=l 

CXD s Iq —1 s 

+ J 2 J 2 '" Y ^bio(7roo)“*+'“^vro6,^j/i'’*+i o5“*+i(P^.(Z))7r(0 

ai = 16i=l flj+i=2 6j+i=l 


(50 S S S 

+ YY'"Y Y ^b,bi+y‘+^ o g{Pp.{l))7r{l) 

ai = l bi=l bj=lbj^i=l 

^ Ti + rs + rs. 


Next we will analyze the boundness of ri,r 2 and Fs one by one. 

OO S OO S 


< EE--- E E h^^^HU)7ril + l) 

ai=l bi=l aj+i = lo 6j+i = l 

= VFi, 


where the inequality is derived from Lemma [2] and VFi is a bounded matrix. 


OO S Iq —1 S 

s EE- E E vrfe,o(7roo)“'+^ o P 0 .{l))Tr{l) 

ai=l bi=l a.j+i=2 bjj^i=l 
oo S Iq —1 S 

= EE--- E E 

ai=lbi=l aj^i=2 bjj^i=l 

oo s Iq—2 s ^z+i“l 

+ E E'" E E E A'Q{A'yw{l) 

(ii=lbi=l <2j+i=2 fej_j_i = l i=0 

oo S Iq —1 5 

+ EE--- E E 7rfe;o(vroo)“*+'"^vrob,+i^'’'+^ 

ai=l bi=l a.j_|_i=2 bjjj.i = l 

= r '2 + + hFs. 
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It is straightforward to verify that W2 and W3 is bounded. 

r, < EE- E 

ai = lfei = l 6 j+i=l 

00 s s 

= EE--- E + 

(11=1 bi = l bj^i=l 

00 S S ^j + 1 — 1 

+ E E'" E ’"«.« E A‘Q{A‘yw{i) 

<11=1 bi=l bj^i=l i=0 

00 s s 

+ E E'" E {lA A'<‘)1J,|.4 «<')]•) (A‘'«)',r(0 

(ii=l bi=l bj^i=l 

^ r'^ + Wi + W 5 , 

where W4 and VI 5 can be readily shown to be bounded. In summary, 

E[P,+i] < r'+ r' + y, forj>l, (14) 

where V = Wy + W2 + VI 3 + W4 + Ws- By a similar argument, we impose (HID to take equality 
and take vectorization. From (fT3|) and (fT4|) . we can calculate E[vec(P/ 3 ,^J] recursively as follows 

E[vec(Pft^J] = r(Pi^)'4'vec(Eo) + T(Pi^)^-i4'vec(0;_i) + • • • + vec(0o). 

where 0o,..., 0z-i are the functions of Q, A, PT^^^’s and are bounded for V is bounded, 


T = \ !,...,! ] ( 8 ) /„ 2 xn 2 


numbers 


and 'k = [V’l,..., 'ip'sl' G with 


lo _ 

= {A® Ay [ ^ vro( 7 roo)“i“ 27 robi(yl“i + P(“i)(:7(“i)) 0 (A'^^ + 

a\=2 


+ 


7 rb,(yl + A'(i)C)®(yl + pWc) , for i G {!,..., s} 


Therefore, E[vec(P^j^J] is bounded as / —>■ 00 if p{Hk) < 1. By some basic algebraic manipula¬ 
tions, one obtains that E||P^;^J| is uniformly bounded if p{Hk) < 1, which completes the proof. 
□ 


The stability condition in Theorem [T] is difficult to test. In the following, we provide an 
equivalent condition. In view of this result. Theorem [1] can be recast as an LMI feasibility 
problem. As for the conversion to LMIs using Schur complements, we refer readers to [20] for 
details. 
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Proposition 1 The following statements are equivalent: 


(i). There exists K = ,..., where ’s are matriees with compatible dimen¬ 

sions, such that p{Hk) < 1, where Hk is defined in (fTTI) ; 

(ii). There exist Xi > 0,..., > 0 and Ki ,..., K\^_i such that 

S lo —1 

{TToof-^A^iA^ + KiC^^'>)XfiA^ + KiC^^^)*{A^y 

i=l 1=2 

s 

+ J2^ij^HA + KiC)Xi{A + KiC)*{A^y < Xj, for all j G 5/{0}.(15) 

i=l 

Proof, (i) {ii) Since p{Hk) < 1, we have 

(/ - Hk)-^ = i + Hk + {Hk)^ + • • • . (16) 

We define a linear space H" as 

VjG5}. 

Then define a norm on as 

s 

II^IU = 5^11^*11 

i=l 

for any H = [Hi,..., Hg] € H”. For any G = [Gi,..., Gs],H = [Hi,..., Hg] G , we say that 
G H (and G >- H) if Gi > Hi (and Gi > Hi) for all i G {1,... , s}. 

Since {I — Hk) is nonsingular and vectorization is a bijective mapping, for any H G 
there exists a unique matrix X = [Xi,..., X^] G such that 

vec{H) = (/ — HK)yec{X). (17) 

In what follows, we shall show X G H". From (1161) . we have 

vec(X) = {I — HK)~^vec{H) 

OO 

= J](77K)Vec(77). 

i=0 

Taking the inverse mapping of vectorization gives X ^ 77, implying X G H”. Similarly, by 
taking the inverse mapping of vectorization on (113, we have 

S lo —1 

H, = Xj-Y, J^(7roo)'-"X^(X^ + X«C«)X,(X' + X«C«)*(X^)' 

i=l 1=2 

s 

- Y + K^^'>C)Xi{A + K^^^G)*{A^y, for all j G 5/{0}. 

i=l 
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where this claim follows as asserted. 


{ii) {i) Dehne an operator Ck — ■ ■ ■, ^K,s) '■ m+ 1HI+ as 

S lo —1 

i=l 1=2 

s 

+ Y + K,C)H,{A + KiCYiA^y,. ( 18 ) 

i=l 

where K = [Ki ,..., iir|^_i], and H = [Hi ,..., Hg] G H”. It is evident that CK{oiH) = aCK{H) 
for any a G M, and that Ck{G) ^ Ck{H) for G H. From the hypothesis of (ii), we conclude 
that there exists a /r G (0, 1) such that Ck{X) :< //X, where X = [Xi,..., X^]. In addition, for 
any given Hq G , there always exists an r > 0 such that Hq A rX. Therefore, for /c G N, 

C%{Ho) A rdkiX) A rfiCYiX) ^ ^ rfi’^X, 

which leads to ||£^(Ffo)||* < As A: —>■ oo, we have limfc_^oo ||£^(Ffo)||* = 0- Note that 

vec {Ck{Hq)) = FfA'vec(Ffo)- Combining all the above observations, we have 

lim {HK)^yec(HQ) = 0, 
fc—>-oo 

which implies p{Hk) < 1- This completes the proof. □ 


4 Comparison with pT 


In this part, we compare our result with those in m and show the advantages of ours. Recall 
that the sufficient condition in m is p($) < 1 where 


$ 4 


d?p+Yi^ooy~"4^^ci 


lo-l 




1=2 


diag(||Af,...,||A^f) , 


with P, Q being defined in (fTT]i and A min^(() ||Xd) + 


4.1 Invariance with Respect to Similarity Transformations 

Theoretically, a state variable transformation (i.e., a similarity transformation from a linear sys¬ 
tem {A, B,C, D) to {S~^AS,S~^B,CS,D) through the nonsingular matrix S does not change 


13 







4.2 Conservativity Comparison 


the stability considered in this work. However, different state variable transformations may gen¬ 
erate opposite conclusions from the stability condition given in HZ!. The invariance of stability 
behavior with respect to state variable transformations can be reflected well from the stability 
conditions presented by this work. 

Proposition 2 Let S € be nonsingular. Suppose there exists K = ..., 

where ’s are matrices with compatible dimensions, such that p{Hk) < 1, where Hk is defined 
in (fTT]l for {A, C). Then, there always exists K = ..., such that < 1, 

where Hj^ is defined for {A,C) = {S~^AS,CS) in accordance with (fTT]l . 

The proof follows from Proposition [Tj and direct calculation. We use the following example to 
illustrate this idea. 


Example 1 Consider the system 

A = 


1.3 0.3 

0 1.2 


C=\l 1 


Q = I 2 x 2 o,nd R = 1, and the bounded Markovian packet-loss process with transition probability 
matrix given by 


n 


0.6 

0.2 

0.2 

0.8 

0.1 

0.1 

0.8 

0.1 

0.1 


(19) 


From flTj , we have d^^^ = 1.2200 and p{^) = 0.7352 < 1. Let 


S = 


1 5 
0 1 


For the system {A,C) = {S ^AS,CS), we have = 1.3632 and /?($) = 1.5202 > 1. 


4.2 Conservativity Comparison 

The stability condition given in this work is less conservative compared with that in m , since 
the latter condition implies the former one. To show this, we need the following proposition. 
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4.2 Conservativity Comparison 


Proposition 3 Define 


^K = 


lo-l 

diP + ^(7roo)'-^d;Q 
1=2 


diag(||A|p,...,||A'in 


where P and Q are defined in (fTT]l and di = |p, and K = [K^^\..., 

with ’s of compatible dimensions. If there exists K such that p{^k) < 1? then p{Hk) < 1- 

Proof. If treating a scalar as the Kronecker product of two other scalars, similar to Proposition[Tl 
we can claim that, if and only if p{^k') < Ij there exists a vector 


X = [xi ,... ,Xs], 


where Xj > 0 for all j € 5/{0}, such that 

S Iq —1 S 

^TTioTTOj '^{TTooy~‘^di\\A^\\‘^Xi + '^TTijdlWA^\\‘^Xi < Xj. 
i=l 1=2 i=l 

The submultiplicativity and subadditivity of a matrix norm result in the following inequality 

S lo —1 

II J^VTioVTo,- Y,{^ooy-^XiA^{A^ + 


2=1 


1=2 


+ ^T:ijXiA\A + K^^^C){A + K^^'^Cy{A^)' < xj, for all j G 5/{0}. (20) 
2=1 

Let Xj = Xjinxn- Then we obtain from (I20p that 

S lo —1 

J^vTiovroj- '^{TTooY-^AyA^ + [A^' 

2=1 1=2 

s 

p'^T^ijAyA + K^^^C)Xi{A + K^^^C)*{Ay’ < Xj. 

2 = 1 

Therefore p{Hx) < 1, which completes the proof. □ 

In virtue of Proposition [31 it is evident that p(<h) < 1 implies p{Hk*) < 1, where K* = 
[Kf,..., with Kf = argminj^(i) 

Example 1 (cont'd) IPe continue to consider Example [I] with an alternative transition proba¬ 
bility matrix 


n 


1 — 


0.6 

0.2 

0.2 

0.6 

0.2 

0.2 

0.6 

0.2 

0.2 
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From Theorem 2 in we obtain p{^) = 1.4704 > 1. By solving an LMI feasibility problem 
using the cvx in Matlab, we see that our Theorem\^ still holds with a group of feasible variables 



0.1081 

0.0243 


-0.8079 

II 

><1 

II 

0.0243 

0.1042 

and K = 

-0.5914 


If we consider the transition probability matrix which only allows the maximum length of 
consecutive packet losses to be 1, i.e., 


n2 


0.6 0.4 

0.8 0.2 


then /?(4>) = 0.49 < 1 and the condition in our Theorem holds. When we increase tth in 112 
from 0.2 to 0.5, one can verify numerically that p($) > 1 while Theorem{l\ of this paper still 
holds. 


5 Conclusion 

We have considered the bounded Markovian packet-loss process model and the notion of the 
peak-covariance stability for the Kalman filtering. A sufficient stability condition with bounded 
Markovian packet losses was established. Different from that of this stability check can 
be recast as an LMI feasibility problem. Then we compared the proposed condition with that 
of |17] . showing that our condition prevails from at least two aspects: i) Our stability condition 
is invariant with respect to similarity state transformations, while the previous result is not; 
2) More importantly, onr condition is proved to be less conservative than the previous one. 
Nnmerical examples were provided to demonstrate the effectiveness onr result compared with 
the literature. 
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